Tight-binding modeling and low-energy behavior 
of the semi-Dirac point 



S. Banerjee\ R. R. P. Singh^, V. Pardoi-^, and W. E. PicketttI 

^Department of Physics, University of California, Davis, CA 95616 and 
^ Departamento de Fisica Aplicada, Universidad de Santiago de Compostela, E-15782 Santiago de Compostela, Spain 

(Dated: June 8, 2009) 

We develop a tight-binding model description of semi-Dirac electronic spectra, with highly 
anisotropic dispersion around point Fermi surfaces, recently discovered in electronic structure cal- 
culations of V02/Ti02 nano-heterostructures. We contrast their spectral properties with the well 
known Dirac points on the honeycomb lattice relevant to graphene layers and the spectra of bands 
touching each other in zero-gap semiconductors. We also consider the lowest order dispersion around 
one of the semi-Dirac points and calculate the resulting electronic energy levels in an external mag- 
netic field. We find that these systems support apparently similar electronic structures but diverse 
low-energy physics. 



The ability to prepare graphene (single graphite 
sheets) [ij has spurred the study of electronic behavior 
of this unique system, in which a pair of Dirac points 
occur at the edge of the Brillouin Zone Bands ex- 
tend linearly (referred to as "massless Dirac") both to 
lower and higher energy from point Fermi surfaces. This 
unusual behavior requires special symmetry and non- 
bonding bands. In bilayer graphene [3|, the linearly dis- 
persing bands become quadratic, while retaining many 
of the symmetry properties. Spin-orbit coupling in these 
systems should lead to a gapped insulator in the bulk. 
Gapless modes remain at the edge of the system, pro- 
tected by topological properties and time-reversal sym- 
metry. This new state of matter, insulating in the bulk 
and metallic at the edges, has been called a topological 
insulator [4, 5]. 

Point Fermi surfaces also arise in gapless semicon- 
ductors in which the bands extend quadratically ("mas- 
sively" ) from a single point separating valence and con- 
duction bands @|. However, these systems are, generi- 
cally, not topological insulators. It has been argued that 
in HgTe quantum wells, where s and p bands overlap 
each other at the F point as a function of well thickness, 
Dirac-like spectra can also arise with exotic topological 
properties. Due to the enhanced spin-orbit coupling in 
these materials, a state of matter exhibiting quantum 
spin Hall effect, has been predicted and observed |8|. 

Recent developments in the synthesis of controlled 
nanostructures, heterojunctions and interfaces of tran- 
sition metal oxides represent one of the most promising 
areas of research in materials physics. While several re- 
cent studies of oxide interfaces have focused on the polar- 
ity discontinuity that can give rise to unexpected states 
between insulating bulk oxides, including conductivity 
0, , magnetism 1 1 , orbital order jT3| , even supercon- 
ductivity [l3|, unanticipated behavior unrelated to po- 
larity can also arise. The V02/Ti02 interface involves 
no polar discontinuity, but only an open-shell charge 
and local magnetic discontinuity, according to the change 
dP across the interface. 



It was recently discovered [14| that a three unit cell 
slab of VO2 confined within insulating Ti02 possesses a 
unique band structure. It shows four symmetry related 
point Fermi surfaces along the (1,1) directions in the 2D 
Brillouin zone, in this respect appearing to be an ana- 
log to graphene. The dispersion away from this point is 
however different and unanticipated: a gap opens linearly 
along the symmetry line, but opens quadratically along 
the perpendicular direction. The descriptive picture is 
that the associated (electron or hole) quasiparticles are 
relativistic along the diagonal with an associated "speed 
of light" vp, as they are in graphene in both directions, 
but they are non-relativistic in the perpendicular direc- 
tion, with an effective mass m. Seemingly the laws of 
physics (energy vs. momentum) are different along the 
two principal axes. The situation is neither conventional 
zero-gap semiconductor-like, nor graphene-like, but has 
in some sense aspects of both. This kind of spectra was 
found to be robust under modest changes in the struc- 
ture. 

Here, we develop a tight-binding model description of 
this semi-Dirac spectra. We find that a three-band model 
is needed, which can be downfolded to two bands at low 
energies. A variant of the model, with only two bands, 
gives rise to anisotropic Dirac spectra, where one has 
linearly dispersing modes around point Fermi surfaces, 
with very different "speed of light" along two perpen- 
dicular axes. A common feature of the various systems 
discussed above are point Fermi-surfaces. From a device 
point of view, for example in thinking of p-n or p-n-p 
junctions, these systems may share common qualitative 
features. The actual dispersion, which would give rise to 
different density of states, may control more quantitative 
differences. However, a more fundamental difference may 
be in their topological properties. [1, [l5, 3, 13] 



We begin with a 3-band tight-binding model of spinless 
fermions [corresponding to the half-metallic VO2 trilayer 
although the system could also be nonmagnetic (spin de- 
generate)], on a square-lattice, defined by the Hamilto- 
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with 63 >> £1,62, so that we have two overlapping 
bands 1 and 2, with no couphng between them. Instead, 
they couple through the third band, by a coupling which 
changes sign under rotation by 90 degrees. Such a cou- 
pling can be shown to arise by symmetry between d and 
s orbitals, for example. The important aspect is that the 
coupling vanishes along the symmetry line, allowing the 
bands to cross (they have different symmetries along the 
(1,1) line). Now, since the third band is far from the 
Fermi energy it can be taken as dispersionless. Further- 
more, without affecting any essential physics, we take 
t\ = —t2 = t and Ai = A2 = t' . Thus, in momentum 
space the Hamiltonian becomes a 3 x 3 matrix: 

^eife Vk\ 
H=\ l2k Vk\ (2) 
Vk 14 £3/ 

where the dispersions and coupling are given by 

£ife = El + 2t(cos kx + cos ky) 
£2fe = £2 - 2t(cos kx + cos ky) 
Vk = 2t' {cos kx — cos ky) 

Using the fact that orbital 3 is distant in energy, the 
three-orbital problem can be downfolded to a renormal- 
ized two orbital problem which becomes (neglecting a 
parallel shift of the two remaining bands) 



H = 




(3) 



The eigenvalues Ek± of H as given by 




£2fe) = 



£3 



(4) 



With some (not very stringent) restrictions on ei — £2 
to ensure that the uncoupled bands actually overlap, the 
two bands touch only at the point ksd along the (1,1) lines 
where £ik — £2k, otherwise the two bands lie on either 
side of the touching point (the Fermi energy). When 
the 2x2 Hamiltonian is expanded around the semi-Dirac 
point ksd it becomes 



H 




VFq2 ql/2m\ 
ql/2m -VFq2) 



(5) 



where 92 and qi denote the distance from ksd along the 
(1,1) symmetry direction, and the orthogonal (1,1), re- 
spectively. The Fermi velocity vp and effective mass m 
can be related explicitly to the tight binding model pa- 
rameters, and also calculated by standard ab initio tech- 
niques. The dispersion relation is that found for the three 
layer slab of VO2 trilaycr in Ti02 at low energy, 



Eq± ^ ±^ {ql/2m)-^ + {vFq2) 



(6) 



For comparison, the graphene dispersion relation is 
= ±VF\/q\ -\- f?!- A plot of the low-energy disper- 
sion of the model giving rise to a semi-Dirac point in 
the 2D Brillouin Zone is shown in Fig. [T] For the VO2 
trilayer,T4] this dispersion holds up to 10-30 meV in the 
valence and conduction bands. 




semiDirac point 
band structure 



Anisotropic Dirac point 
band structure 



FIG. 1: On the left, the plot shows the low energy band 
eigenvalues Eq± in a region near Ef for the semi-Dirac point. 
On the right is the same plot for the anisotropic Dirac point. 



A few observations can be made at this point. First, 
if the original bands 1 and 2 were simply coupled by the 
same anisotropic mixing Vk (without any third band in 
the picture), then anisotropic Dirac points (rather than 
semi-Dirac points) occur along the (1,1) directions. This 
dispersion is also shown in Fig. [TJ This type of two- 
band situation should not be particularly unusual, hence 
Dirac points in 2D systems are probably not as unusual 
as supposed, i.e. they are not restricted to graphene nor 
are they restricted to high symmetry points. 

While the constant energy surfaces of our model may 
appear to be elliptical (the common situation; the Dirac 
point has circular FSs), they are actually quite distinct. 
As E — >0 the velocity is constant in one direction and 
is y2rnE in the other; the FSs vanish as needles with 
their long axis perpendicular to the (1,1) direction. This 
can be seen in Fig. [2l showing the constant energy sur- 
faces for electron doping according to our model, show- 
ing the 4 semi-Dirac points in the tetragonal k^j-kj, Bril- 
louin zone. The density of states (DOS) n(E), which 
is constant for effective mass systems and goes as \E\ 
for graphene, is proportional to y^jS] at a semi-Dirac 
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point. When doped, the density of carriers will follow 
n{Ep) cx \Ef\^/'^ behavior. 




FIG. 2: Plot showing the Fermi surfaces for electron doping 
that derive from the low energy excitation spectra of the semi- 
Dirac point at the F — A4 direction in the k^-kj, plane of the 
square Brillouin zone. 

Another observation is that the same bands Eq± can 
be obtained from related but distinct low-energy models, 
such as 

and 

Although the bands resulting from H2 and H3 are the 
same, the eigenfunctions are different and are intrinsi- 
cally complex for H2 and H3 unlike for the specific semi- 
Dirac point we discuss. 

One of the issues of most interest to such systems is the 
behavior in a magnetic field. Making the usual substitu- 
tion q p+ I A with momentum operator p and vector 

potential A, we find the Landau gauge A = B{—X2,0, 0) 
to be the most convenient here. First, however, we note 
that the characteristics of the two directions, the mass 
m and velocity vp, introduce a natural unit of momen- 
tum po — mvp and length Xo — fi/Po, and of energy 
Co = ■mvp/2. Introducing the atomic unit of magnetic 
field Bo such that fisBo = 1 Ha, and the dimension- 
less field b = B/Bo, units can be scaled away from the 
Hamiltonian by defining for each coordinate x 

X2 = (^)'/^Xoi2, (9) 

and similarly for xi. Here 7 is the dimensionless ratio 
of the two natural energy scales: 7 — iiBBo/{mvp/2). 



Under this scaling 

Pi + -Ai =pi~ -Bx2 ^ Po{lbf'\pi - ^2) (10) 
c c 

where pi,xi are conjugate dimensionless variables, etc. 
Thus all possible semi-Dirac points (all possible m and vp 
combinations) scale to a single unique semi-Dirac pointy 
with the materials parameters determining only the over- 
all energy scale. There is no limiting case in which the 
semi-Dirac point becomes either a Dirac point or a con- 
ventional effective mass zero-gap semiconductor. For the 
case of trilayer VO2, 7 does not differ greatly from unity 
0- 

Shifting X2 io u = X2 — Pi, with conjugate dimension- 
less momentum p, the Hamiltonian in a field becomes 

H - 2e, {jbf/'' [p a, + ^u^ a^] (11) 

= 2eo (76)'/' h. 

The energy scale is much larger than for conventional 
orbits though smaller than in graphene so the VO2 
trilayer may display an integer quantum Hall effect at 
elevated temperature as does graphene [isj . 

A scalar equation for the eigenvalues can be obtained 
from /i^. Introducing the operator Q = p + iu'^/2, the 
eigenvalues of are Q^Q and QQ^ , giving the mathe- 
matical problem 

1 

Q^QMu) = i--r^ + -u''- u)Mu) = ^iMu)- (12) 
au^ 4 

The equation for QQ^ has the opposite sign of the lin- 
ear term, with identical eigenvalues and eigenfunctions 
related by inversion. Note that every eigenfunction of h 
is also an eigenfunction of , and that although the po- 
tential is negative in the interval (0,4^/^), the eigenvalues 
£^ must be non-negative. 

We have obtained the eigenvalues both by precise nu- 
merical solution and by WKB approximation, finding 
that the latter is an excellent approximation. Initially 
neglecting the linear term in the potential, the WKB 
condition [3] 

.V2e| I — ^ 

J_^,^el--uidu={n+-)7r (13) 

can be solved to give the WKB eigenvalues for the quartic 
potential as 

= [3y|f§y]'^'(- +l^'= l-3765(- +1)^- (14) 

The linear perturbation corrects the eigenvalues only to 
second order, which is significant only for the ground 
state (~0.74 versus the numerical solution of 0.59). The 
WKB error is less than 0.01 for the first excited state 
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FIG. 3: Potential energy function for the one-dimensional 
Schrodinger equation and the resulting quantized energy lev- 
els of h^. The lowest numerical vales of the three energy 
eigenvalues = -t-x/el" are provided. 



metallic ferromagnetic system where time-reversal sym- 
metry is broken. Applications of the VO2 trilayer and 
related semi-Dirac point systems may provide unusual 
spintronics characteristics and applications. 

In conclusion, we have developed a tight-binding model 
description of the semi-Dirac and anisotropic Dirac spec- 
tra relevant to V02-Ti02 multi-layer systems. Our tight 
binding model contains nothing unconventional, indicat- 
ing that semi-Dirac and anisotropic point systems are 
not as rare as has been assumed. The low energy char- 
acteristics of the semi-Dirac point are intermediate be- 
tween those of zero-gap (massive) semiconductors and 
Dirac (massless) point systems. The study of such ox- 
ide nano-heterostructures has only just begun and they 
clearly promise a number of diverse electronic structures 
and novel phases of matter. 

This project was supported by DOE grant DE-FG02- 
04ER46111 and by the Predictive Capabihty for Strongly 
Correlated Systems team of the Computational Materials 
Science Network. V.P. acknowledges financial support 
from Xunta de Galicia (Human Resources Program). 



and gets successively smaller for higher eigenvalues. We 
see then that the semi-Dirac system has eigenvalues in 
a magnetic field which scale as 5^/'^ and increase as 
(n + ^)^^'^ as n gets large. Both aspects lie between 
the behaviors for conventional Landau levels (linear in 
B, proportional the n+ ^) and the Dirac point behavior 
(proportional to VB n), as might have been anticipated. 
Some low-lying eigenvalues of are shown in Fig. [3] 
against the potential well. Note that there is no zero- 
energy solution as in the graphene problem. 

Another way in which Dirac spectra can arise on a 
square-lattice can be motivated in terms of the model of 
Bernevig et al. for HgTe quantum wells. In their 
model the two bands crossing each other have s and 
p characters respectively. Thus the interband hopping 
term changes sign under reflection. This can lead to 
a (sin kx ~t~ isin/uy) coupling between the bands. Note 
that in this model, only a single Dirac point can occur 
and it must be at /c = 0, when the two bands touch 
each other at that point. In contrast, in the models dis- 
cussed here, there are four symmetry related semi-Dirac 
(or anisotropic Dirac) points whose location can vary con- 
tinuously along the symmetry axis (1,1), with changes in 
band parameters. A feature unique (so far) to the VO2 
trilayer system is that point Fermi surface arises in a half 



[1 
[2: 

[s: 
[4: 

[s: 
[6: 

[7 
[8 
[9 

[10: 
[11 

[12 

[13 
[14 

[15 
[16 
[17 

[18 
[19 



Electronic address: I wepickett ©ucdavis . edu | 
K.S. Novoselov et al, Nature 438, 197 (2005). 

A. H. Castro Neto et al, Rev. Mod. Phys. 81, 109 (2009). 
T. Ohta et al, Science 313, 5789 (2006). 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005) . 

L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
I. M. Tsidilkovski, Gapless Semiconductors (Akademie- 
Verlag Berlin, Berlin, 1988). 

B. A. Bernevig et al, Science 314, 1757 (2006). 
M. Konig et al, Science 318, 766 (2007). 

A. Ohtomo et al. Nature 419, 378 (2002). 

A. Ohtomo and H.Y. Hwang, Nature 423, 427 (2004). 

A. Brinkman et al, Nat. Mater. 6, 493 (2007). 

R. Pentcheva and W.E. Pickett, Phys. Rev. B 74, 035112 

(2006) . 

N. Reyren et al.. Science 317, 5842 (2007). 

V. Pardo and W. E. Pickett, Phys. Rev. Lett. 102, 

166803 (2009). 

D. J. Thouless et al, Phys. Rev. Lett. 49, 405 (1982). 
F. D. M. Haldane, Phys. Rev. Lett. 93, 206602 (2004). 
J. E. Moore and L. Balents, Phys. Rev. B 75, 121306 

(2007) . 

K. S. Novoselov et al. Science 315, 1379 (2007). 

C. M. Bender and S. A. Orszag, Advanced Mathemati- 
cal Methods for Scientists and Engineers (McGraw Hill, 
1978), p. 523. 



